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Abstract. We present results from a calculation of the transition temperature in 
QCD with two light (up, down) and one heavier (strange) quark mass as well as for 
psj ! QCD with three degenerate quark masses. Furthermore, we discuss first results from 

an ongoing calculation of the QCD equation of state with almost realistic light and 
strange quark masses. 
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1. Introduction 

Lattice calculations provide a unique first principle approach to the thermodynamics 
of QCD. They have the potential to produce definitive answers for basic bulk 
thermodynamic properties as well as the phase diagram of QCD. To reach this ambitious 



goal one has to control lattice discretization errors that lead to systematic deviations 



of lattice calculations from continuum physics, as well as finite volume effects that that 
can obscure the thermodynamic limit. 

We present here results from a detailed analysis of the transition temperature in 
QCD with an almost physical value of the strange quark mass and several values of 
the light quark mass. This allows to perform an extrapolation to the chiral limit 
as well as to the physical point. Calculations have been performed for two different 
values of the lattice cut-off which allows to estimate the lattice discretization errors 
and perform an extrapolation to the continuum limit. The flavor dependence of these 
results are examined through an analogous analysis of QCD with three degenerate quark 
flavors. We also report on results from an ongoing study of the QCD equation of state 
with almost physical light and strange quark masses. All these calculations have been 
performed with an C(a 2 ) tree-level improved gauge action and an improved staggered 
fermion action (p4fat3) that also improves bulk thermodynamic observables at 0(a 2 ). 
For further details on the calculational set-up we refer to [U E]. 

2. The transition temperature 



In order to determine the transition temperature and its dependence on the light quark 
mass we have performed calculations with a fixed value of the bare strange quark mass 
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Figure 1. The transition temperature (left) calculated on lattices with temporal 
extent N T = 4, 6 and the static quark potential (right) in units of r calculated on 
16 3 x 32 lattices. V(r) has been calculated for a wide range of gauge coupling (3 in 
which the lattice cut-off varies by a factor 4. This range also covers the the relevant 
parameter range of the finite temperature calculations. 



and several (degenerate) values for the two light quark masses that correspond to light 
pseudo-scalar masses (m ps ) ranging from m ps ~ 150 MeV to m ps ~ 500 MeV. For 
non-zero values of the quark masses the transition to the high temperature phase of 
QCD is signaled by peaks in response functions. We have analyzed Polyakov loop and 
chiral susceptibilities to locate the transition point for the different quark mass values 
on lattices with temporal extent N T = 4 and 6 pQ. The spatial extent of the lattice 
has been varied between N a /N T = 2 and 4. We generally observe only a weak volume 
dependence of the transition point and also find that the above susceptibilities as well 
as quark number susceptibilities in the light and strange quark sector [3] yield, within 
statistical errors, consistent results for the transition point. 

In order to convert the results for the transition temperature to physical values 
one needs to set the scale through the calculation of zero temperature observables. 
Unlike hadronic observables in the light quark sector, which require a careful handling 
of chiral extrapolations, are known to be difficult to determine and in some cases are 
influenced by large discretization errors, the parameters of the static quark potential 
have been found to show only a mild dependence of the lattice cut-off p£J. We thus have 
performed extensive studies of the static quark potential, V(r), on zero temperature 
lattices and extracted the scale parameter r defined through the slope of V(r), i.e. 
(r 2 dV(r)/dr ) = 1.65. Using the results obtained for r in lattice units at the 

\ / r=r'o 

crossover couplings, which have been determined in finite temperature simulations on 
lattices with temporal extent N T = 4 and 6, the transition temperatures can then be 
calculated as T c r = r /N T a. Results for T c r and the static quark potential in units of 
ro are shown in Fig. [TJ Results for T c ro have then been extrapolated to the continuum 
limit using an ansatz that incorporates the quark mass dependence of T c through the 
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Figure 2. The trace anomaly, (e — 3p)/T 4 , calculated in (2+l)-flavor QCD on lattices 
with temporal extent N T = 4 and 6, respectively. Fig.2(right) shows the resulting 
pressure from Eq. [3] as well as earlier results obtained for 3- flavor QCD (see text) . 

calculated values of the lightest pseudo-scalar mass as well as a quadratic dependence 
on the lattice cut-off, 

+ A{m ps r Q ) d + B/N 2 T . (1) 

Here d has been varied between d = 1 and 2 to allow for uncertainties in the extrapolation 
to the chiral limit flj. The continuum extrapolation is shown as the lower curve in 
Fig. deleft). For the transition temperature at the physical point we find, 

T c r = 0.457(7)^ <S> T c = 192(7) (4) MeV , (2) 

where we have used r = 0.469(7) fm [5] to set the scale for T c . 

A similar analysis has been performed for QCD with three degenerate quark 
masses [2]. In this case the transition becomes first order for small but non-zero 
values of the quark mass. We find for the transition temperature in the chiral limit 
T c r = 0.419(9) while in the same limit the transition temperature in (2+l)-flavor QCD 
is T c tq = 0.444(6) (12). We note that T c ro in 3-flavor QCD is about 5% smaller than in 
(2+l)-flavor QCD. Moreover, the difference between the transition temperature in the 
chiral limit and at the physical point is only about 3%. This is a direct consequence of 
the weak quark mass dependence of the transition temperature, which has been observed 
earlier [6]. A similarly weak quark mass dependence has also been observed in studies 
performed with another version of an improved staggered fermion action (asqtad) on 
lattices with temporal extent N T = 4, 6 and 8 [7j. This analysis, however, lead to a 
10% smaller transition temperature, T c r = 0.402(29) than the value reported here. To 
resolve the origin of this discrepancy needs further investigations. 

3. The QCD equation of state 

We currently study the equation of state in (2+1) -flavor QCD [8] on lattices with 
temporal extent N T = 4 and 6 using the same calculational set-up used for the study of 
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the transition temperature. Calculations are performed along a line of constant physics 
characterized by a fixed strange pseudo-scalar mass m js ~ 660 MeV and a fixed ratio 
of TUss and the kaon mass, m^/m^ — 1.3. The light pseudo-scalar (pion) mass on this 
trajectory is about 220 MeV. 

In Fig. [2] we show results for the trace anomaly, 0^(T) = e — 3p in units of T 4 and 
the pressure (p) which is deduced from 0^(T) using standard thermodynamic relations, 

We also included in Fig. fright) results for the pressure obtained earlier [9] for 3-flavor 
QCD with larger quark masses that have been held fixed in units of T. For these earlier 
calculations the light pseudo-scalar mass was about 770 MeV in the transition region 
and would have been about 2 GeV in the region around T ~ 3T C ! 

The results presented here for (e — 3p)/T A and pjT A in (2+l)-flavor QCD are 
entirely consistent with results obtained with the asqtad action [10] . In particular, 
results for p/T 4 obtained with the p4fat3 and asqtad actions, which both have a 0(a 2 ) 
improved ideal gas limit, show little cut-off dependence in the entire temperature regime 
analyzed. This is in contrast to calculations performed with a variant of standard 
staggered fermions [11] which unlike p4fat3 and asqtad actions lead to large cut-off 
distortions of the high-T limit on lattices with temporal extent N T = 4 and 6. We also 
note that we find no evidence for a significant broadening of the transition region when 
going closer to the continuum limit. Whether this happens for larger N T , as a recent 
analysis of the transition temperature on lattices with temporal extent up to N T = 10 
may suggest [12], has to be examined in more detail. 
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